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> ■ Abstract 

a^ : 

CN . I present measurements of the projected rotational velocities of a sample of 

I 100 early B-type main-sequence stars in the Large Magellanic Cloud. This is the 

• first extragalactic study of the distribution of stellar rotational velocities. The 

. sample is drawn from two sources: a sample derived from the vicinity of the 

main-sequence turnoff of young clusters (ages l-3xl0^yrs), and a sample from 
the general field. I find the cluster population exhibits significantly more rapid 
rotation than that seen in the field. I have drawn analogous Galactic cluster and 
O I field samples from the literature. Comparison of these samples reveals the same 

^ I effect. I propose the observed difference between cluster and field populations can 

d ■ be explained by a scenario of evolutionary enhancement of the surface angular 

^ , momentum over the main-sequence lifetime. A comparison is made between the 

cluster and field populations of the LMC and the Galaxy in order to explore the 



o 

Oh' 



X 



^ • effects of metallicity. This shows that the LMC stars are more rapid rotators 



than their Galactic counterparts. 
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1 Introduction 

This paper aims to establish the distribution of projected rotational velocities (vsini) 
for a sample of early B-type stars in the Large Magellanic Cloud (LMC). 1 explore the 
effects of the parameters age and metallicity on the observed distribution of rotational 
velocities. The motivation for this study arises from recent stellar evolutionary models 
that incorporate stellar rotation. This modeling shows that the rotational history is 
vital to the interpretation of the evolutionary phase and elemental abundances seen in 
massive stars. Heger, Langer & Woosley (2000) and Meynet & Maeder (2000) show 
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that axial rotation can bring about in a natural way the degree of extension to the 
convective core in massive stars as discerned from cluster H-R diagrams (e.g. Keller et 
al. 2000, 2001a) and the mass-luminosity relation of Cepheids (Keller & Wood 2002). 

Models incorporating rotation also predict mixing of internally processed material 
to the surface during main-sequence (MS) evolution (the faster the initial rotation the 
greater the degree of mixing). Standard evolution models (e.g. Bressan et al. 1993), 
on the other hand, do not. In the Magellanic Clouds (MCs), nitrogen has been shown 
to be over-abundant in many supergiants of spectral types B to K (e.g. Venn 1999; 
Lennon et al. 1996) and numerous MS B stars (Rolleston ct al. 1996; Korn et al. 2000; 
Korn et al. 2002) which is indicative of mixing of CNO cycled material. The dispersion 
in N abundance is significantly larger than that seen in the Galaxy (Mc Erlean et al. 
1999; Venn 1995; Gies & Lambert 1992). This has been tentatively attributed to the 
presence of a greater number of rapid rotators within the lower metallicity environs of 
the MCs. 

At present, direct observational studies of the dependence of stellar rotation on 
metallicity are limited. Within the Galaxy the work of Burki & Maeder (1977) exam- 
ined the variation of mean vsini of early B stars with galactocentric distance, but the 
range of the latter was only several kpc which, given the galactic metallicity gradient, 
is insufficient to draw firm conclusions regarding the role of metallicity. Recent work by 
Royer et al. (2003) has revisited this work extending it to a broader sample of Galactic 
clusters with inconclusive results. 

Indirect observational evidence for more rapid rotation of massive stars in metal- 
poor sites hke the MCs is derived from examination of the population of Be stars. Be 
stars show Balmcr emission arising from a circumstellar disk of material surrounding 
a rapidly rotating stellar photosphere. They form a population with which to trace 
stellar rotation (assuming that the onset of the Be phenomenon itself is not governed 
by metallicity but rather the underlying distribution of rotational velocities). Maeder 
et al. (1999) examine the fraction of Be stars (i.e. N(Be)/N(B)-|-N(Be) ) against metal- 
licity for a sample of clusters with relevant ages (1-3x10^ years, corresponding to the 
maximum occurrence of Be stars). Maeder et al. find a strong anti-correlation between 
the Be fraction and metallicity. This leads the authors to conclude that rotation rates 
are higher amongst metal-poor stars. 

The role of age in determining the resultant distribution of vsini was first suggested 
by Abt & Hunter (1962) in their study of the brighter members of the Trapezium and 
Pleiades clusters. This was elaborated upon by Wolff et al. (1982) and Guthrie (1984) 
in their studies which contrasted the distribution of rotational velocities seen in young 
clusters, older associations and field stars. The distribution of rotational velocities 
in the older associations and the field share a high proportion of very slow rotators 
compared with the young cluster sample. 

Keller et al. (1999, 2000, 2001b) have used the Be star population to explore the 
age dependence of the distribution of rotational velocities. We have observed that 
the Be fraction amongst the young cluster population rises towards the luminosity of 
the MS terminus, whereas in the field it retains a relatively uniform distribution in 
luminosity. This difference between cluster and field populations can, we argue (Keller 
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et al. 2001b), be brought about from an evolutionary enhancement in the rotational 
velocity over the later portion of the MS lifetime. 

It is the aim of this paper to present direct measurements of the projected rota- 
tional velocities for a sample of early B-type MS stars within young cluster and field 
environments of the LMC, together with analogous Galactic data from the literature, 
to investigate directly the dependence of vsini on both age and metallicity. Section 2 
presents our observations and method for the determination of projected rotational ve- 
locities. Section 3.1 compares the distribution of rotational velocities in the cluster and 
field environments of the Galaxy and the LMC. Section 3.2 compares the Galactic and 
LMC cluster populations to explore metallicity effects. Finally in Section 4 1 present a 
discussion of the physical causation of the observed dependence of the distribution of 
projected rotational velocities on age and metallicity. 

2 Observations and data reduction 

I have defined two sample sets: (1.) A cluster sample - drawn from the population in 
the vicinity of three young populous clusters in the LMC, namely, NGC 1818, NGC 
2004 and NGC 2100. Our targets were selected from the WFPC2-based photometry 
of Keller et al. (2000). On the basis of this photometry I have selected those stars on 
the MS with 14.0 <V < 16.0. (2.) A field sample - from the vicinity of NGC 2004 
and 1818. These regions were selected to form the basis of our sample because of their 
richness and apparently small differential reddening (Keller, Wood & Bessell 1998). 

Medium resolution (0.6A/px) spectra of the sample were obtained using the Double- 
Beam Spectrograph (DBS) on the SSO 2.3m telescope between 4-6 March 1999, 25-16 
November 1999 and 10-12 January 2001. The spectra consist of two simultaneously 
recorded non-overlapping segments: blue (3800-4800A) and red (6200-6800A). The 
observations made use of aperture plates machined to locate the target stars. In general 
the density of the targets precluded the use of slits for each object, rather holes of 
appropriate size for typical seeing at the site (2" ) were made. Sky correction was latter 
made through the use of shts at the extremities of the field. Wavelength cahbration 
was made through CuAr spectra interleaved with the observations. 

2.1 Determination of rotational velocities 

To determine the projected rotational velocity, vsini, I have used a minimisation 
technique which locates an optimal match between the observed spectrum and rota- 
tionally broadened synthetic spectra. The grid of synthetic profiles of stellar H Balmer 
and Hel lines by Gonzales Delgado & Leitherer (1999) is used. This grid of models 
spans 50000> T^s >4000 K and 0.0> logg > 5.0. For T^s > 25000 K NLTE models 
are used to compute the synthetic spectra, for cooler stars, Kurucz LTE models are 
used. 

To create the rotationally broadened spectra the synthetic spectra were first con- 
volved with the instrumental response function and then with a broadening function 
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which is based upon that of Gray (1976): model stellar surface is divided into 40 seg- 
ments and the emergent doppler-shifted intensity from each is calculated. This is then 
integrated to yield the flux at each wavelength. 

This approach assumes that the resultant line profile is independent of viewing 
angle. Physically, the surface gravity (and hence temperature) of a rotating star varies 
with latitude as a result of gravity darkening (von Zeipel 1924). As a consequence the 
contribution of the rapidly rotating cooler equatorial regions to the observed spectrum 
is reduced. Townsend et al. (2004) present a survey of the behaviour of hue width with 
vsini including viewing angle dependent gravity darkening. For instance, a star rotating 
at 95% of the critical velocity (i.e. the velocity at which the equatorial escape velocity 
drops to zero) seen pole on using the method we apply here leads to an underestimation 
of VeSini by 20%. The underestimation of VeSini is a strong function of Ve such that 
by Ue=300kms~^ for a B2 star this has fallen to a few percent. The neglect of gravity 
darkening underlies all rotational velocity standards and previous studies of rotational 
velocity distributions. For this reason the present work maintains this simplification. 

I have used a wavelength range of 3990A to 4550A for the purposes of finding 
an optimal match. The spectra is fit as one contiguous piece. An optimal solution is 
sought in vsini, Tes, logg parameter space. The spectra of a series of rotational velocity 
standards (Slettebak et al. 1975) of early to mid B spectral type were obtained. Figure 
1 shows the comparison between the vsini determined by means of our parameter 
fit procedure and that of Slettebak et al. (ibid). Figure 2 shows a typical fit to the 
rotational velocity standards. Figure 1 does not show any systematic difference between 
these results and those of Slettebak et al. 

The signal-to-noisc ratio (S/N) of the spectra of our LMC sample is low, a typical 
3hr exposure results in a S/N=30 on a V=15.5 star. On the basis of simulations 
(discussed below) 1 set a lower limit to the S/N of our sample of 20 {V ~ 16). Table 1 
& 2 reports the vsini for our sample. Some typical fits are shown in figure 3. 

The uncertainties in our estimated vsini are a function of both S/N and vsini. I 
have carried out a series of trials using model spectra which are rotationally broadened 
and then degraded to a given S/N. 1 have then examined the vsini returned by our 
routine from repeated trials. Slow rotators arc perturbed by a low S/N. A S/N of 
20 introduces an uncertainty of ±50 kms^^ at a vsini of kms^^ and systematically 
returns a higher vsini. This overestimation of vsini at low S/N results primarily from 
the radial velocity correction which must be made for each star during fitting. At low 
S/N the line-centers of the significant H and He lines are less distinct and an optimal 
fit may occur with cither a slightly higher or lower radial velocity than the underlying 
systemic radial velocity. A S/N degraded line seen off-center appears as a broader line 
and this results in the systematic shift of the slowest rotators to higher vsini. 

As vsini is increased the importance of this effect diminishes. At vsini=100 kms~^ 
and S/N=20 the uncertainty in vsini falls to ±25 kms~"^ with no significant systematic 
shift. The uncertainty grows again to ±75 kms^-*^ at vsini=350 kms~^ (S/N=20) as the 
line depth is reduced but again with no significant systematic shift. This behaviour 
is summarised in figure 4. By describing in detail the systematic behaviour of our 
technique at low S/N it is possible to statistically remove it as described in section 3.2. 
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Figure 1: Comparison between projected rotational velocities of Slettebak et al. (1975, 
open circles) and those of the present study. Solid points show the positions of B15 
and B30 in NGC 2004 from Korn et al. 2000. 
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Figure 2: Examples of the fits between artificially broadened theoretical profiles (thick 
line) and observed stellar spectra for rotational velocity standards; a. ( Cru (B2.5V, 
t'sim=70kms~^) and b. 6 Cru (B2IV, vsim=140kms~^) 
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Figure 3: Examples of rotationally broadened line profile fits for LMC field stars. Plot 
a. shows fp2_8 (^=13. 83 and t;sim=190kms~^) and b. shows fp3_9 (1^=15.89 and 
i'sim=50kms~^). 
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Table 1: vsini for the sample of main-sequence stars drawn from the LMC clusters 
NGC 1818, 2004 and 2100. V magnitudes are taken from Robertson (1974). vsini 
with additional colon indicate those objects for which determination of vsini have 
considerable uncertainty (±75 kms~^). 
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Table 2: vsini for the sample of main-sequence B stars drawn from the LMC field. V 
photometry is that of the authors. 
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Figure 4: Histograms of the vsini returned from our fitting routine from three models 
with different input vsini (0, fOO, 350 kms~^ resp.) and S/N. The thick hne shows the 
distribution recovered from a S/N=50 and the thinner hne from S/N=20. 

3 Distribution of rotational velocities 

3.1 Comparison of cluster and field populations 

Figure 5 compares the distribution of rotational velocities from LMC cluster and field 
samples. The non-parametric two-sample Kolmogorov-Smirnov (K-S) test is conducted 

to test the validity of the proposition that both samples are derived from the same 
parent distribution. Such a test reveals that the cluster and field are unlikely to arise 
from a similar parent distribution at 2a significance (Prob.(D>D"'")=4%; D"'"=0.25, 
Ni=49, N2=51). 

To extend our analysis to a discussion of the role of metallicity in the distribution 
of rotational velocity I have constructed an analogous sample from the Galactic field. 
Our LMC sample has a limit of V ~16 imposed by the attainable S/N. Due to the high 
inclination angle and the essentially disk-like morphology of the LMC we can regard 
the sample as effectively volume-hmited. Assuming a distance modulus to the LMC of 
18.45 our sample extends to a My—S. This corresponds to a spectral type of B2 on 
the MS (Zorec and Briot 1991). 

The Galactic field sample is drawn from the Bright Star Catalogue (Hoffeit & 
Jaschele 1982:BSC) and the Supplement to the Bright Star Catalogue (Hoffeit et al. 
1988:SBSC). I have drawn from the catalogue stars of spectral types B0-B2 and of 
luminosity class HI-V. The BSC-I-SBSC sample is magnitude-hmited at V—7.1. In 
order to compare with the volume-limited sample from the LMC we first must reduce 
the apparent number of stars within each spectral type and luminosity class division 
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Figure 5: The distribution of rotational velocities from the 49 LMC cluster stars (solid 
line) and 51 field stars (dotted) shown on the left as a histogram and to the right as 
the non-parametric cumulative distribution. 

to a volume-limited distribution. This must be done otherwise the most luminous 
stars (i.e. early spectral type and high luminosity class) will receive a disproportionate 
weighting (Mamquist bias). 

I shall consider a common volume defined by the BOIII stars. The sample-filling 
factor e for each spectral type and luminosity class is given in Zorec and Briot (1997, 
their table 4) for the BSC-I-SBSC sample. The corrected histogram distribution of 
vsini, hcorivsini) is constructed from the sum of the histograms for the individual 
spectral types and luminosity classes, C, scaled by e: 

hcor{vsmi) = ■ ^i^i '■ vsini) (1) 

i 

where e ranges from 4.5 for B2V to 1.0 for BOIII. The net result is a distribution 
which is strongly weighted by the input of the intrinsically most numerous members, 
that is, those of B2V. 

A distribution of projected rotational velocities was drawn from a sample of young 
Galactic clusters (see table 3). The chosen clusters are rich clusters with analogous 
ages to those of the LMC cluster sample. Stars within 2 magnitudes of the MS turnoff 
(B0-2V-III) were selected for the distribution. 

The distributions of rotational velocities for the two samples are compared in figure 
6. The results of a K-S test on the above distributions shows that the Galactic cluster 
stars are significantly faster rotators than their field counterparts (Prob.(D>D+)=2%, 
D+=0.165, Ni=310, N2=151). 
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Table 3: Input to the Galactic cluster vsini distribution. 
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Figure 6: The distribution of rotational velocities from the Galactic cluster (solid 
line) and field (dotted) shown on the left as a histogram and to the right as the non- 
parametric cumulative distribution. 
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Figure 7: The distribution of projected rotational velocities for the Galactic cluster 
sample (dotted line) and for the LMC cluster sample (solid line). The dashed line 
is the distribution for the Galactic cluster sample after simulation of the systematic 
effects of low S/N (see text). 

3.2 Comparison of Galactic and LMC populations 

Let us now compare the distributions of vsini in both the LMG and Galactic samples 
to investigate the possible role of metallicity in determining vsini. Figure 7 shows the 
cumulative distribution functions of the cluster and Galactic field samples (solid and 
dotted lines respectively). As discussed in section 2.1, the accuracy of our determina- 
tions of vsini are critically limited by the obtainable S/N. To statistically account for 
the systematic effects imposed by low S/N I have taken the spectral type, vsini data 
of the Galactic field sample and produced a set of synthetic spectra with a distribu- 
tion of S/N mimicking that of the LMC sample and have then determined the vsini 
distribution of the population as for the LMC sample. This is shown in figure 7 as the 
dashed line. As expected, the net effect is to shift a proportion of the slowest rotators 
to higher vsini. 

The resulting K-S test between the LMC cluster sample and the S/N degraded 
Galactic cluster sample rates the probability of these two samples been alike at 6% 
(D+=0.22, Ni=49, N2=151). This is suggests that the average rotation rate is higher 
amongst the LMC cluster stars. Table 4 summarises our results. 

4 Discussion — the role of evolution and metallicity 

Our observations present us with two important conclusions: 1. Young cluster stars 
(BO-2 V-III within 2 mag. of the MS terminus) are faster rotators than a similar 
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Table 4: mean vsini for the four samples discussed above. 



Location 



mean vsini (kms ^) 



Galactic field 
Galactic young clusters 
LMC field 



85 
116 
112 
146 



LMC young clusters 



luminosity range in the field, and 2. LMC cluster stars are faster rotators than Galactic 
cluster stars within the same luminosity range at just under a 2a significance. 

To account, firstly, for the observed difference in rotational velocity between the 
cluster and field populations 1 propose a scenario of evolutionary rotational enhance- 
ment. This scenario utilizes the one clear difference between the cluster and field 
populations, that is, the age spread within each population. 

Consider a spin-up phase which occurs over an interval of the MS lifetime towards 
the end of the MS. The small age spread in the cluster population at the luminosity 
of the cluster MS terminus places the majority of these stars within the interval of 
spin-up. The field population, on the other hand, possesses a more uniform spread of 
ages at a given luminosity, from the ZAMS to the MS terminus (this can be confidently 
said for such rapidly evolving stars) . Hence in the field we see the time-average of the 
MS rotational velocity. 

The evolution of rotational velocity has been examined by Endal & Sofia (1979) 
and more recently, by Hcgcr, Langcr & Woosley (2000) and Meynet & Maeder (2000). 
These models explicitly follow the radial exchange of angular momentum during the 
course of stellar evolution. 

The studies of Hegcr, Langer & Woosley (2000) and Meynet & Maeder (2000) focus 
on stars of >20Mq which is more massive than that considered in our sample (12-5Mq; 
Bressan et al. 1993). A major distinction between massive stars and those considered 
in our sample is the dominance of mass loss on the course of MS evolution. Their high 
mass-loss rates lead to the removal of large amounts of angular momentum. Mass-loss 
is enhanced amongst fast rotators that approach the critical angular velocity . The 
removal of large amounts of mass from a >2OM0 star results in a decline in the angular 
velocity over the course of MS life. 

The study of Endal & Sofia (1979) focuses on the rotational velocity evolution of 
a 5M0 star. Here, the mass-loss rates are much reduced. The authors find that stars 
commencing their MS lives with relatively slow angular velocities remain slow rotators 
throughout the course of the MS evolution (although there is a marked increase in 
angular velocity of all stars during the core contraction phase at the exhaustion of 
central hydrogen burning this phase is far too short to account for any observable 
increase in the number of rapid rotators). By contrast, stars commencing their lives 
with an angular velocity greater than 60-80% of the critical angular velocity spin up 
towards the critical velocity over a moderate fraction of the MS lifetime. 
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The largest masses considered here are of the order of 12M0. Such stars have a 
mass loss rates 1.5 dex lower than those of a 20Mo star (de Jager et al. 1988) and 
consequently the MS evolution is not dominated by mass-loss. Hence we can assume 
the behaviour of the present sample is best described by the models of Endal & Sofia. 

In a cluster, therefore, those stars closer to the MS turnoff are more likely to be 
rapid rotators since a proportion of them will have been spun up by evolution. Less 
luminous stars that have not evolved as far through the MS phase will not have entered 
the spin up phase. This process could provide the mechanism required to explain the 
more rapidly rotating sample seen in the cluster population in the vicinity of the MS 
turnoff. 

Our second finding, that the LMC cluster and field samples exhibit more rapid 
rotation than their Galactic counterparts indicates a metallicity effect. Our sample 
is insufficient to determine whether the metallicity dependence resides in the initial 
distribution of rotational velocities or if the magnitude of the evolutionary spin-up 
discussed above increases in lower metallicity environments. Our observational sample 
is necessarily limited by the small sample size attainable with our instrumentation. To 
verify the findings presented here we await the extension to a larger LMC (and SMC) 
sample. This is currently underway utilising the FLAMES multi-object spectrograph 
on the VLT. 

5 Summary 

In this paper 1 have presented a comparison of the rotational velocities of a sample of 
stars within the field and from a selection of young clusters in both the Galaxy and, 
from our own observations, the LMC. I have presented the distribution of 100 LMC 
early B-type main-sequence stars. This represents the first extragalactic study of the 
distribution of stellar rotational velocities. 

It is found that the early B-star population of young clusters (1-3x10'' yrs.) in 
both the Galaxy and the LMC exhibits more rapid rotation than the field population. 
I propose this can be explained by a scenario of evolutionary enhancement of the 
surface angular momentum brought about by angular momentum redistribution over 
the main-sequence lifetime. 

A comparison of a sample of field and cluster stars drawn from both LMC and 
Galactic environments shows that the LMC stars a more rapid rotators than their 
Galactic counterparts. The origin of this metallicity dependence is as yet unknown. 
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